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We present the first 3-dimensional high resolution hydro- dynamical simula- 
tions of the interaction between the hot ionised intra-cluster medium and 
the cold interstellar medium of spiral galaxies. Ram pressure and turbu- 

nous galaxies like the Milky Way within 100 million years. These mechanisms 



featureless disky galaxies contain no atomic gas and show no signs of recent star- formation 



naturally account for the morphology of SO galaxies, the rapid truncation of 
star formation implied by spectroscopic observations, as well as a host of ob- 
servational data on the HI morphology of galaxies in clusters. 

CO 

Crucial observational evidence for the hierarchical formation of structure in the 
universe is the dramatic evolution of galactic morphologies in dense environments over 
the past 5 billion years ([UH). The key puzzle that remains to be solved is the origin 
of the large population of lenticular (SO) galaxies found in nearby clusters These 

Or 
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The Hubble Space Telescope revolutionised our view of the universe by revealing that 
distant galaxies appeared different from the local population. In contrast to local clusters, 
^ ■ high resolution imaging of distant clusters led to the spectacular finding that young clusters 

of galaxies are filled with spiral galaxies (||J4],||) and contain almost no lenticular (SO) 
galaxies, whereas the ratio of luminous ellipticals to lenticulars (SO) increases by a factor 
of five between a redshift z=0.5 and the present-day (||). SO's can be characterised by their 
thick featureless disks that show no evidence for recent star-formation and the increase in 
their population appears to be countered by a similar decrease in the number of luminous 
late-type spirals in clusters. The data suggest that a transformation between these galaxy 
types is taking place as a direct consequence of the cluster environment. 

Three mechanisms have been proposed that can lead to morphological transforma- 
tion between galaxy classes. Mergers will transform disks to spheroidals (|6|,[/]), but are 
only effective in low density environments (^). Gravitational tidal interactions between 
cluster galaxies can naturally account for the observed evolution of the faint end of the 
luminosity function and the transformation of small disks to faint spheroidals (^). How- 
ever, more massive bulge dominated systems are stable to tidal disruption (|io|). Although 
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the resulting disk thickening from tidal heating suppresses spiral features and causes a 
morphological similarity to SO's, neither of these two processes suppress star-formation. 

In addition to disk thickening, a mechanism that actively extinguishes star-formation 
is crucial since the stellar populations of SO galaxies are old and their spectra indicate that 



star-formation was abruptly halted several billion years ago (11.12,13). A slow decline in 
the star formation rate, such as expected from the exhaustion of a reservoir of cold gas, is 
unable to explain the strongly enhanced Hydrogen absorption lines seen in many distant 
cluster galaxies. 

A candidate mechanism was proposed by Gunn & Gott (|l4|) over two decades ago. 
Their simple force-balance estimates suggested that the motion of galaxies through the hot 
ionised intra-cluster medium (ICM) creates a "ram-pressure" that could potentially strip 
away significant amounts of gas from disks. However, a full description of this mechanism 
must include complex turbulent and viscous stripping (|l^) at the interface of the cold and 
hot gaseous components as well as the formation of bow-shocks in the ICM ahead of the 
galaxy. Although these processes have been cited over 1000 times in the literature, their 
effectiveness and efficiencies have received little theoretical investigation. 

The ram pressure is proportional to p ICM Vg al , therefore the infalling galaxy suffers 
most gas loss at its pericentric passage, where its velocity can be as large as 3000 km s _1 
and the intra-cluster gas density approaches 3000 /i^atoms m~ 3 . Gunn Sz Gott's (14) 



order of magnitude estimates suggest that the typical gas disk would be stripped down to 
~ 5 kpc, a radius confirmed by smoothed particle hydro-dynamic simulations that follow 
just the ram pressure process (|i~6|). However, this leaves 50% of the original HI confined 
in the disk which would continue to form stars for several billion years. Although this is 
a significant reduction in the star formation rate, it does not explain the absence of any 
recent star formation in cluster SO galaxies. 

We have performed the first high resolution three dimensional numerical simulations 
of these hydro-dynamical processes to accurately address the efficiency of stripping and the 
timescale on which it occurs. Previous work has been either in 2D with spherical galaxy 
models ( p^Jl^ ), or in 3D using a code that could not model viscosity and turbulence (16). 



Our parallel computer code (19) uses high-resolution shock- capturing techniques to follow 
the fluid dynamic equations allowing us to observe the full complexity of the ram-pressure 
and turbulent/viscous stripping processes. We can also follow the shocks that penetrate 
the interstellar medium (ISM) and ICM and the thermo-dynamical evolution of the ICM 
and stripped galactic gas. Because it is fully three dimensional, our code is not forced to 



preserve the cylindrical symmetry of the galaxy (|1 

We construct a self-consistent equilibrium model galaxy with stellar disk and bulge 
components designed to resemble a luminous spiral similar to the Milky Way or An- 
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dromeda. The stars are embedded within a dark matter halo constructed such that the 
total rotational velocity of the disk is a constant 220 km s" 1 (20). The real ISM is complex 
gaseous medium formed by a cold diffuse HI component and dense molecular clouds (MC) 
with temperatures Tmc ~ 10 2 K. The diffuse gaseous disk is constructed by specifying 
the density, and velocity at each grid cell according with (|20|), and constant temperature 
Tism = W A K- the lower threshold in our simulations. The MC's are typically three order 
of magnitude denser than the HI component and have sizes of the order of several parsecs. 
Even though our numerical code is highly optimised and is running on state of the art 
parallel hardware, the maximum resolution that we can achieve is ~ 100 parsecs, therefore 
we cannot resolve individual MC's in our simulations. However, MC's are so small and 
dense that they will remain unaffected by the stripping processes (pl|). Initially, the ICM is 
considered as an uniform medium with constant density and temperature Tjqm = 10 8 K. 
Both the ICM and ISM are treated as ideal fluids with adiabatic exponent 7 = 5/3. Their 
evolution is described by the hydrodynamic equations that are integrated using the nu- 
merical techniques described above. The stellar and dark matter components are evolved 
using an N-body Particle-Mesh code. All components are coupled gravitationally through 
Poisson's equations which is solved using a fast 3D FFT method. No cooling has been 
considered since the timescales are so short (p2[). At our best (standard) resolution we use 
512 3 (256 3 ) cells across a cubic region centered on the galaxy with length 64 kpc, thus our 
nominal resolution is ~ 100(200) parsecs. 

We have simulated different infall geometries, velocities and ICM gas densities and 



ISM structure (23). Two of these simulations are pictured in Figure 1, which show a 
galaxy moving face on and nearly edge on through the core of a rich cluster like Coma 
- p ICM = 2.6 x 10 3 h^Q 2 atoms m -3 - at a velocity of 2000 km s —1 . We plot only the 
ICM and ISM components of the simulation to highlight the gas dynamical processes - 
the stellar disk remains unperturbed. The outer gas disk of the infalling galaxy is rapidly 
stripped away and forms trailing streams of warm HI in pressure equilibrium with the 
ICM. Turbulence and viscous stripping ablate the gas disk even further (24). 

We find a rich structure of shocks. The most obvious is the prominent curved bow 
shock that propagates through the ICM ahead of the galaxy, heating the ICM from a 
temperature of 10 8 to 5 x 10 8 degrees ( p5"l , p6| ). A complex series of cross shocks occur in 
the rarefied hot medium behind the galaxy, which may be visible as a wake of enhanced X- 
ray emission (p5|). Although not apparent on this scale, a second shock is driven through 
the ISM of the galaxy, raising the internal pressure by over two orders of magnitude. 
The efficiency and time scale of star formation are strong functions of the ambient ISM 
pressure (|27|). As the shock propagates through the galaxy, it may promote the collapse 
of molecular clouds, briefly enhancing the galaxy's star formation rate. 

Whereas previous work treated the ISM as a smooth disk of HI, in reality, the ISM 
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has complex multi-phase structure, filled with bubbles, shells and holes ranging in size from 
a few parsecs to a kpc fl27|). Furthermore, the inner couple of kpc of most bright spiral 
galaxies are extremely deficient in HI. The nearest and best studied galaxy is Andromeda 
which has over 100 HI holes and a central region of radius 2 kpc devoid of neutral gas 
(|28|). Our simulations show that this structure makes the disk much more susceptible to 
viscous stripping. As the ICM streams through these holes, it ablates their edges and 
prevents stripped gas falling back on to the centre of the galaxy. This is an important 
difference between previous work that claimed significant replenishment from stripped 
material. When we model our disk on Andromeda's we find that these processes lead to 
the removal of the entire diffuse HI component on a timescale of 100 million years. Even 
if the HI holes contain a large quantity of molecular hydrogen locked within MC's, due to 
their small covering factor they would not affect the removal of the HI component of the 
gaseous disk (21). Combined with galaxy harassment, this process has all the properties 
required to explain the rapid transformation between spiral and SO galaxies seen in distant 
clusters. 

It is important to stress that the effectiveness of ram-pressure stripping does not 
depend on the galaxy moving face-on through the ICM. We find that galaxies inclined 
twenty degrees to the direction of motion suffer as much stripping as face on encounters. 
Since the orientation of the disk remains fixed as the galaxy orbits through the cluster, a 
galaxy moving edge on at first pericentric passage is likely to be inclined to the direction 
of motion at some other point of its orbit (^) . Only a small fraction of galaxies will orbit 
with their disks aligned exactly with the orbital plane. 

The timescale for stripping is very short compared to the orbital timescale, therefore 
galaxies will only rarely be observed at the moment of stripping. Evidence for ongoing 
ICM/ISM interaction may be found by observing compression of the leading contours of 
disk HI emission, and wings or tails of HI behind the galaxy (3^,31). In Figure 3 (see also 
32) we compare two recent and striking examples of these processes with snapshots of our 
simulations, showing that the observed morphologies are well reproduced. After passage 
through the cluster, disks will be either HI deficient compared to field spirals, or they 
may lack neutral hydrogen altogether. Galaxies in this latter class would subsequently be 
identified as SO's. HI maps of galaxies in the Virgo and Coma clusters show few galaxies 
with any HI, although some galaxies in the outer parts of the cluster are extremely HI 



deficient (p4|. 33 J35|) . The efficiency of stripping explains why such galaxies are so rare. 

Our simulations demonstrate the importance and effectiveness of the ram-pressure 
and transport processes. The interaction between the ICM and ISM removes the entire 
HI component as well as any diffuse reservoir of hot gas within its dark matter halo. 
But what happens to the molecular clouds in the infalling galaxy? These clouds are so 
dense and small that they cannot be removed by the ram-pressure of the ICM and are 
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not resolved in our simulation. To understand the fate of the molecular gas we need to 
consider the cycle of star formation in a little more detail. In a quiescent field galaxy, 
molecular clouds are continually being disrupted by the star clusters formed within them 
and they are subsequently reformed from the ambient diffuse HI. Within clusters, galaxies 
are swept clear of HI and this cycle is broken: the disruption of the clouds is not balanced 
by the condensation of new self-shielding molecular complexes. Current models suggest 
the lifetime of large molecular clouds is less than a few tens of millions of years 
would therefore expect the ram pressure to lead to a decline in star formation on the same 
timescale as the disruption of the molecular clouds. This scenario predicts that cluster SO 
galaxies will not contain molecular gas (|37|). 

Furthermore, the pressure increase in the ISM may create a burst of star-formation, 
as the counter shock compresses existing molecular clouds within the galaxy (|38|). This 
would explain the prevalence of strongly enhanced Hydrogen absorption lines in many 
distant cluster SO galaxies, and is supported by the spectacular H Q emission seen from 
several nearby galaxies as they are stripped ( |39l , f40| ). 

This mechanism leads to a population of post-starburst SO galaxies, preferentially 
populating the central regions of rich clusters. However, ram-pressure is not the only effect 
at work in clusters of galaxies. Longer timescale processes, such as stellar evolution, tidal 
forces and 'harassment' tend to thicken the stellar disk and enhance the relative impor- 
tance of the galaxy's bulge. Taken together, these processes can explain the remarkable 
transformation of cluster galaxies and the dramatic evolution in the galaxy populations 
of dense environments. SO's in the field most likely form via a minor merger or the ac- 
cretion of a satellite (^). This formation mechanism is markedly different from what we 
are proposing in clusters, and suggests that there will be observable photometric, spectral 
and kinematical differences between field and cluster SO's. 
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Ps(R,z) = J^2 z exp(-R/R s )sech 2 (z/z s ), 

where R s = 3.5 kpc, z s = 0.35 kpc, M s = 5.6 x 10 10 M . (iv) Gas disk, 

Pg( R , z ) = exp(-R/R g )sech 2 (z/z g ), 

being R g = 3.5 kpc, z g = 0.35 kpc, M g = 1.4 x 10 10 M . The total masses of the 
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r coma r coma \r coma 

10 3 fill) 2 atoms m~ 3 ) . The ICM velocities used were 1000 kms^ 1 and 2000/cms~ 1 . 
The orientation of the galaxies moving through the wind were varied from face-on 
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FIGURES CAPTIONS 



Figure 1 The evolution of the gaseous disk of a spiral galaxy moving face on (left 
column) and inclined 20 degrees to the direction of motion (right column) through a 
diffuse hot intra-cluster medium. Each snapshot shows the density of gas (5 = p/p ICM ) 
within a 0.2 kpc slice through the center of the galaxy and each frame is 64 kpc on a side. 
Note how rapidly the disk material is removed - within 100 million years 100% of the HI is 
lost. We do not show the stellar disk, bulge or dark matter halo which remain unaffected 
by the loss of the gaseous component. The box size is 64 kpc and the hydro grid has 256 3 
cells. 

Figure 2 Mass loss as a function of time for the model including ten little holes and 
inhomogeneous density. We plot the evolution of the gaseous mass within a cylindrical 
slice of 25 kpc radius and thickness 2 kpc centred on the center of mass of the stellar disk. 
Initially, ram-pressure stripping dominates the gas loss process and the entire outer disk is 
removed in a very rapid timescale. Viscous and turbulent stripping operates continuously, 
but over a longer timescale, resulting in a roughly linear rate of mass loss. We stop this 
simulation after 130 million years by which time 97% of the gas disk has been removed. 

Figure 3 Observational evidence (left panels) for ram-pressure processes compared with 
our hydro-dynamical simulations at different epochs (right panels). The first panel shows 
a HI map of NGC 7421 which shows wings of gas being pushed back by its motion through 
a diffuse ionised medium ([H]). The second panel shows the HI deficient galaxy NGC 4548 
orbiting in the Virgo cluster ([30|). The remaining gas has a ring- like morphology very 
similar to our simulation after 50 million years. 



Next Figures would be shown in http:/ /www. scienceonline.org as an extra 



material 



Figure 4 Observational evidence (left panels) of trails of gas compared with our sim- 
ulations at different epochs (right panels). The first panel shows the radio continuum 
(contours) brightness distribution superposed over a gray scale representation of H inten- 



sity in galaxy 97073 in cluster A 1367 (see Figure la in 39). The second panel shows 
spectacular radio jets that have been swept backwards for tens of kpc by motion of the 
galaxy through the ICM of a rich cluster at z=0.3 (radio data courtesy of R. Ivison, AW. 
Blain and I. Smail). Note that the radio emission may be unrelated to the stripping pro- 
cess, but the morphological appearance is very similar to the trails of stripped gas in our 
simulations. 
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